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Abstract 

The S 2 valued wave map flow on a Lorentzian domain RxE, where S is any flat 
two-torus, is studied. The Cauchy problem with initial data tangent to the moduli space 
of holomorphic maps £ — > S 2 is considered, in the limit of small initial velocity. It is 
proved that wave maps, in this limit, converge in a precise sense to geodesies in the 
moduli space of holomorphic maps, with respect to the L 2 metric. This establishes, in a 
rigorous setting, a long-standing informal conjecture of Ward. 

1 Introduction 

Wave maps are the analogue of harmonic maps in the case where the domain is Lorentzian. 
They satisfy a semilinear wave equation which has been heavily studied, in the simplest non- 
trivial case of S 2 target space, as a model PDE system involving a manifold-valued field 
|18j . The wave map equation is particularly interesting in the case where the domain is 
(R x Y,,dt 2 — <?s), with (S,ps) an oriented Riemannian two-manifold. In this case, the static 
wave map problem is conformally invariant, so static solutions on £ = R 2 in particular have 
no prefered scale: they can be dilated to any size without changing their energy. This suggests 
that time- dependent solutions with initial data close to a static solution might collapse and 
form singularities in finite time, an issue which has been heavily studied both numerically 
and analytically mainly for E = R 2 , N = S 2 , within a certain rotational equivariance class 
[U [10J LT2J [16] . It turns out that finite-time collapse can indeed occur, and the collapse mecha- 
nism is stable (to small perturbations of the initial data) at least within the equivariance class 



This paper addresses a different analytic issue from singularity formation, namely the 
validity of the geodesic approximation to wave map flow. To motivate this, one should think 
of wave maps R x £ — > S 2 as (formal) critical points of the action functional 

S=jdt{T-E), where T = ±J |^| 2 , E= l -jjd^\ 2 . (1.1) 
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It follows (from Noether's Theorem) that they conserve the total energy T + E. A rather 
general argument of Lichnerowicz [11] shows that for any map : E — > S 2 of topological 
degree nGZ (subject to suitable boundary conditions, if S is noncompact), E > 47r|n|, with 
equality if and only is is ± holomorphic. So holomorphic maps, if they exist, minimize 
potential energy in their homotopy class. Let us denote by M n the moduli space of degree n 
holomorphic maps S — > S 2 . Consider a wave map <f)(t) with 0(0) G M n and 0t(O) G T </i ( )M n 
with ||0t(O)||i2 small. By conservation of E + j\\(pt{t)\\ 2 L2 , one expects that 0(t) will stay 
close to M n , on which E attains its minimum value, for as long as the solution persists. This 
led Ward to suggest [2E], in the specific case S = IR 2 , that such wave maps should be well 
approximated by the dynamical system with action S, but with </>(£) constrained to M n for all 
time. Since E is constant on M n , this constrained system is equivalent to geodesic motion on 
M n with respect to the L 2 metric (obtained by restricting the quadratic form T to TM n ). A 
similar approximation had previously been proposed by Manton [T3] for low energy monopole 
dynamics, and the geodesic approximation is now a standard technique in the study of the 
dynamics of topological solitons [TH] . 

Geodesic motion on M2 (for E = M 2 ) was studied in detail in [9]. There is a technical 
problem: the L 2 metric is only well-defined on the leaves of a foliation of M n and one must 
impose by hand that <f>(t) remains on a single leaf. This turns out to be ill-justified (it precludes 
singularity formation for n — 1, for example, in contradiction of [16]). This technical deficiency 
is removed if we choose X to be a compact Riemann surface. Here geodesic motion in M n is 
globally well-defined, if incomplete [17j, and the L 2 geometry of M n is quite well understood, 
at least for some choices of £ and n [T3l [191 EH IZE] . 

The question remains: is geodesic motion in M n really a good approximation to wave map 
flow in the adiabatic (low velocity) limit? The purpose of this paper is to prove that it is, 
for times of order (initial velocity) -1 at least in the case where £ is any flat two-torus. More 
precisely, we will prove: 

Theorem 1.1 (Main Theorem). Let M„ denote the moduli space of degree n>2 holomorphic 
maps from a flat two-torus X to S 2 . For fixed O £ M n and (pi G T^ Q M n consider the one 
parameter family of initial value problems for the wave map equation with 0(0) = 0o, 0i(O) = 
e0i, parametrized by e > 0. There exist constants r* > and e* > 0, depending only on the 
initial data, such that for all e G (0,£*] ; the problem has a unique solution for t G [0,r*/£:]. 
Furthermore, the time re-scaled solution 

£ : [0, r„] x S -»■ S 2 , £ (r, p) = <j>(r/e, p) 

converges uniformly in C l to ip : [0, r*] x S — > S 2 , the geodesic in M n with the same initial 
data, as e — >■ 0. 

To prove this we will adapt the perturbation method devised by Stuart to prove validity 
of the geodesic approximation in the critically coupled abelian Higgs and Yang-Mills-Higgs 
models [221 [23] ■ The wave map problem has a key similarity with these gauge-theoretic prob- 
lems, namely a moduli space of static solutions which minimize energy in their homotopy 
class and satisfy a system of first order "Bogomol'nyi" equations. (For wave maps, the Bogo- 
mol'nyi equation is the condition that be ± holomorphic, i.e. the Cauchy- Riemann equation.) 
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Roughly, the idea is to decompose the solution <p{t) as = if>(t) + e 2 Y(t) where ip{i) G M„, 
and control the growth of a suitable Sobolev norm of the error Y(t) uniformly in e by means 
of energy estimates. One concurrently shows that the projected trajectory if>(t) converges to 
a geodesic in M„. 

In comparison with Stuart's work on vortices and monopoles, the situation we study is 
simpler in two respects: we work on a compact domain £ (rather than R 2 or R 3 ), and our 
system has no gauge symmetry. On the other hand, the wave map problem introduces two 
new challenges for the method. 

First, our field is manifold- valued, so it is not clear a priori what the decomposition (f>(t) = 
ifj(t) + e 2 Y(t) really means. In preliminary work on this problem, it was suggested that the 
correct formulation was (j)(t) = exp^n\ e 2 Y(t), where exp : TS 2 —> S 2 is the exponential map 
[5]. In fact, this turns out not to have the analytic properties required by Stuart's method 
(except for rotationally equivariant wave maps). In this paper we isometrically embed S 2 
in R 3 and use the ambient linear structure to project as usual, <fi(t) = ip{t) + e 2 Y{t). This 
choice is simple, but has significant repercussions: Y is no longer tangent to the map ip (not 
a section of ip^TS 2 ), and must satisfy a nonlinear pointwise constraint to ensure that <fi is S 2 
valued. The evolution of Y is governed by a nonlinear wave equation whose (spatial) linear 
part is the Jacobi operator for the harmonic map if> : £ — > S 2 . It turns out that is not 
self-adjoint when acting on non-tangent sections (such as Y). Since self-adjointness of (the 
analogue of) is crucial for Stuart's method, we must devise a way round this: we replace Jf 
by an "improved" Jacobi operator L^, which coincides with on tangent sections, but is self- 
adjoint on all sections, and introduce compensating nonlinear terms into the wave equation 
for Y using the pointwise constraint. Further difficulties result: L^, unlike Jf, does not define 
a coercive quadratic form on the L 2 orthogonal complement of T^M n . We must work instead 
with a weaker near-coercivity property, which turns out to suffice for our purposes. 

Second, while £ is compact, the moduli space M n is not. Of course, the vortex and 
monopole moduli spaces, dealt with by Stuart, are also noncompact, but in those cases, 
moving to infinity corresponds to (clusters of) solitons separating off and escaping to infinite 
separation, a well-controlled process. For wave maps, by contrast, approaching the boundary 
of M n at infinity corresponds to one or more lumps collapsing and "bubbling off". In this 
process, if> becomes singular and both geodesic motion and wave map flow become badly 
behaved. To handle this, we must keep careful track of the position (of if) G M n ) dependence 
of our various estimates, and modify Stuart's a priori energy bound so that we simultaneously 
control the error Y(t), the deviation of ip(t) from the corresponding geodesic, and the distance 
of if>(t) from dooM n . 

It is interesting to speculate to what extent Theorem 1 1.1 1 can be generalized. It is clear that 
the proof presented here generalizes quite easily to the case of a general compact Riemann 
surface, provided n is sufficiently large compared with the genus of S. The reason for restricting 
to the case £ = T 2 is mainly one of presentation: the existence of global cartesian coordinates 
makes it straightforward to define the various function spaces, for example. Generalizing the 
target space is not so straightforward. The wave map flow R x S — > N has the appropriate 
"Bogomol'nyi" form for Stuart's method to apply whenever S, N are both compact kahler 
manifolds (in fact, it suffices for E to be co- kahler). The choice N = CP k , k > 2, is of 
some interest in mathematical physics, for example. But here the reliance on an isometric 
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embedding N cW becomes very problematic. It seems likely that some variant of Theorem 
11.11 does remain true for general compact kahler targets, but proving it would require a rather 
different approach, perhaps along the lines sketched in [24] . 

One should note that Theorem 11.11 gives no information about singularity formation for 
wave maps on R x E because, although there certainly are geodesies ip(r) which hit <9ooM n 
in time r < oo, and the corresponding wave maps do converge uniformly to ip(r) on some 
interval [0,r*], there is no reason to expect t* = r . In fact Rodnianski and Sterbenz have 
shown that singularity formation of equivariant wave maps on R 2 deviates significantly from 
the dynamics predicted by the (suitably regulated) geodesic approximation [16J. Since blow 
up is a (spatially) local phenomenon, these results presumably apply in some form on the 
torus, which would imply r* < r . Nontheless, the geodesic approximation (on compact E or, 
regulated, on £ = R 2 ) predicted finite time blow-up of wave maps in (2 + 1) dimensions, and 
this prediction turned out to be correct. The geodesic approximation also makes predictions 
about the genericity of blow up. It is not hard to prove, for example, that generic geodesies 
on Mi for £ = S* 2 do not hit the boundary at infinity. It would be interesting to see whether 
the full wave map flow has this property (i.e. generic Cauchy data tangent to Mi have global 
smooth solutions). The analogue of Theorem 11.11 for £ = S 2 could provide a starting point 
for proving such results. 

The rest of the paper is structured as follows. In section [2] the moduli space M n of holo- 
morphic maps is introduced and its key property, Proposition 12.11 (existence of a smooth local 
parametrization about any point), established. In section [3], the projection of the wave map 
flow to M n is defined, and the coupled system satisfied by ip and the error section Y is derived, 
equation (13.151) . In section [U some standard functional analytic definitions and results are 
introduced. Our aim here, and in the remainder of the paper, is to make the proof accessible to 
a wide mathematical physics audience, not just experts in PDE. In section a local existence 
and uniqueness theorem for the coupled system (13. 15j) governing (ip, Y) is proved, Theorem 
15.11 Of course, local existence and uniqueness of wave maps in this context is not new; the 
extra, and new, information we obtain here is local existence and uniqueness of the projection 
to M n . This is the engine underlying Stuart's method, and we go through the argument in 
some detail, not only because there are certain new aspects we have to deal with which Stuart 
did not (e.g. preservation of the pointwise constraint on Y), but also because the requirements 
of Picard's method for this proof determine our choice of function spaces, a point which is 
not obvious (to the non-analyst) in Stuart's original applications of the method [221 123]. In 
section [6] the key near-coercivity property of the quadratic form associated with the improved 
Jacobi operator is proved, Theorem 16.51 (roughly, that {L^Y, L^L^Y) L 2 controls the H 3 norm 
of Y). In section [7] energy estimates are established which bound the growth of Y(t). Finally, 
in section [SX the coercivity properties and energy estimates are combined to prove long time 
existence of the solution (i/j,Y), and establish convergence to the corresponding geodesic in 
M n . An appendix presents the proofs of some basic analytic properties of the nonlinear terms 
in the coupled system for (if), Y) which are used repeatedly in section 
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2 The moduli space of static n-lumps 



Static wave maps are harmonic maps : E — > S 2 , that is, solutions of the harmonic map 
equation 

0^ + ra + (|0x| 2 + |0j 2 )0 = O (2.1) 

or, equivalently, critical points of the Dirichlet energy 

E(<P) = \ [ I0,| 2 + W 2 - (2-2) 
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Maps E — > S 2 fall into disjoint homotopy classes labelled by their degree n G Z, which we 
may assume, without loss of generality, is non-negative. An argument of Lichnerowicz [3J p39] 
shows that, in the degree n class, E(<p) > 47m, with equality if and only if is holomorphic. 
Furthermore, all harmonic maps E — > S 2 of degree n > 2 are holomorphic [3J. So the moduli 
space of interest, M n , is the space of degree n holomorphic maps E — > S 2 . Such maps are 
called "n-lumps" by analogy with the case E = C, where the Dirichlet energy density typically 
exhibits n distinct local maxima, which may loosely be thought of as smoothed out particles, 
or lumps of energy. 

The global topology of the space M n is quite complicated, for example, M 2 = [E x 
PSL(2, C)]/(Z2 x Z2), [20]. For our purposes local information will suffice, however. Given a 
smooth variation S of = O G M n we have dE((p s )/ds = at s = (since is harmonic) 
and 

= (V, J,V) = jv ■ J,V. (2.3) 
where V = d s <p s \ s= o G r(0 _1 TS' 2 ) is the infinitesimal generator of the variation and 

J*V = ~V XX - Vyy ~ (|0,| 2 + \<j) y \ 2 )V - 2{<t> x ■ V x + <j) y • Vy)<f> (2.4) 

is the Jacobi operator at the map (ft [25, pl55]. This operator is self-adjoint and elliptic, and its 
spectrum determines the stability properties of <fi. By the Lichnerowicz argument, <fr minimizes 
E in its homotopy class, so spec is non-negative. Given a variation <p s through harmonic 
maps, that is, a curve in M n through (f> = O , its infinitesimal generator V = d s (p s \ s=0 satisfies 
J^V = 0. Hence T^M n C ker J^. For a general harmonic map : M — > N between Riemannian 
manifolds, the converse may be false, that is, there may be sections V G ker C T((f)~ l TN) 
which are not tangent to any variation of through harmonic maps, and in this case the space 
of harmonic maps M — )■ iV may not be a smooth manifold around 0. It is important for us to 
rule out this kind of bad behaviour in our case. More precisely, we need: 

Proposition 2.1. Given any 0o G M n; n > 2, there exists an open set U C M 4n and a smooth 
map ip : U x E — > S 2 such that, 

(i) for each q G U , ip(q, •) G M n , 

(ii) there exists q G U such that 4> = i[)(q , ■), and 

(Hi) ip^ = dip/dq^, [A — 1, 2, . . . , 4n span ker JW g> .) . 
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Proof. Choose any p G S 2 such that both p and — p are regular values of 0o (such p exists by 
Sard's Theorem). Then : X — > S 2 is in M n if and only if s p o0, its image under stereographic 
projection from p, is meromorphic, of degree n, that is, a degree n elliptic function. The most 
general degree n elliptic function is [8] 

( x\( \ - ai ) ■ ■ ■ a{z - a n ) 

(s p o ){z) = X— ( — j— 2.5 

where a is the Weierstrass sigma function, A, a±, . . . , a n , bi, . . . , b n are complex constants, A ^ 0, 
Yl a i = S k mod A and {a^} PI {bj} = 0. Hence, we may parametrize a general point G M n 
by 4n real numbers g M , for example, the real and imaginary parts of X,a±, . . . , 6„_i having 
set b n — ai + • • ■ + a n — b± — ■ • • — 6 n -i- Further, manifestly depends smoothly on q and 
z. Hence we have a smooth map ip '■ U x S — > S 2 satisfying properties (i) and (ii). By our 
choice of p, s p o O has n distinct zeroes and n distinct poles, so {"0 M } at q = g are linearly 
independent sections of i[>(q, •) -1 T5' 2 , and hence, by smoothness, also linearly independent on 
some neighbourhood U of go in U. As explained previously, ip^ span a subspace of ker J^( q ,.), 
so it remains to show that ker has dimension 4n for any G M n . 

It is known [251 pl74] that ker is isomorphic, as a complex vector space, to H°(E,L), 
the space of holomorphic sections of the line bundle L = 0~ 1 T'S' 2 , where T'S 2 denotes the 
holomorphic tangent bundle of S 2 . Since has degree n and T'S 2 has degree 2, L has degree 
2n. Now, by the Riemann-Roch formula [B] 

dim#°(£,L) - dim ^(S, L) = degL = In (2.6) 

since S has genus 1. But, by Serre duality, if 1 (S,L) = H°(H,K ® L*)* where is the 
canonical bundle of E. Now is trivial, so K ® L* has degree — 2n, and hence has no 
holomorphic sections, whence H l (Y,,L) = 0. It follows that ker has real dimension 4n, as 
required. □ 

We can regard q^ as local coordinates on M n . Given the initial data 0i G T^Mn of interest, 
we choose and fix such a ip : U x £ — > S 12 and denote by go £ U and gi G M 4n those vectors such 
that 0o — i/>(<lo, ") an d 0i = gf'0 M (go5 •)■ Here, as henceforth, we use the Einstein summation 
convention on repeated indices. We also choose and fix a compact neighbourhood K of go in 
U. In a slight abuse of notation, we will also use the symbol ip to denote the associated map 
U D K —>• M n C (5* 2 ) E , so -0(g) will denote the holomorphic map z !->■ ip(q,z). We will also 
use ^(t) as shorthand for tp(q(t)), meaning a general curve in M n . 

There is a natural Riemannian metric on M n , the L 2 metric, whose components in the 
local coordinate system g M are 

w«) = = (") 

The associated Christoffel symbol is 

Gl(q) = r a (Tp a ,TpXu) (2.8) 

where 7^ is the inverse metric and ip^ u = d 2 ip / dq^dq u . This is the metric whose geodesies 
approximate wave maps in the adiabatic limit. 
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3 Projection of wave map flow and the coupled system 



The wave map equation for : K x X — > S 2 C M 3 is 

4>tt - 4>xx - 4> yy + (\4>t\ 2 - \4>x\ 2 - \4> y \ 2 )4> = o. (3.1) 
Given e > 0, a small parameter, we decompose as 

= ^ + £ 2 F (3.2) 

where, at each fixed time, ip(t, •) : X — > S 2 is a degree n harmonic map, and Y : E — > R 3 . 
We may think of ip(t) as a curve in M n , the moduli space of degree n harmonic maps, and Y 
as the "error" incurred by projecting (p(t) to ip(t). It is useful to think of Y as a section of 
i/'~ 1 K 3 _, where R?_ = S 2 x R 3 is the trivial R 3 bundle over S 2 , and ■0~ 1 M 3 _ is its pullback to S. 
With this is mind, we refer to Y as the "error section" , and to any Z : S — > R 3 with Z ■ ip = 
pointwise as a "tangent section" (in bundle language, Z is a section of ip^TS 2 C ■^~ 1 R 3 ). 
Clearly, K is not a tangent section (unless K = 0). Since both if) and (f> are S 12 valued, Y must 
satisfy the pointwise constraint 

i\)-Y = - l -e 2 \Y\ 2 . (3.3) 

For a given curve if)(t), if is a wave map then Y must satisfy the PDE obtained by substituting 
f l3~2|) into ([31]), 

Ytt + J^Y = k + ej (3.4) 

where is the Jacobi operator associated with the harmonic map if){t) : E — >■ S* 2 and the 
terms on the right hand side are 

k = -(if) TT + \i(; T \ 2 if)) 

j = 2(W- Y t )if) + e{(\Y t \ 2 - \Y X \ 2 -\Y y \ 2 )iP + (\iP T \ 2 -2ip x -Y x -2iPyY y )Y} 

+2e 2 {i, T ■ Y t )Y + e 3 {\Y t \ 2 ~ l^| 2 - ^y^- (3-5) 

We have here introduced the slow time variable r = et as a book-keeping device. The precise 
expression for j is not important. What matters is its qualitative form: it depends only on 
if>,if) T and Y and its first derivatives, and the dependence is smooth (polynomial, in fact). 

Superficially (\3A\i looks exactly analogous to the corresponding equation in Stuart's anal- 
ysis of vortex dynamics [22], but this is deceptive. As already noted, the Jacobi operator 
is a self-adjoint, elliptic, second order linear operator : T^ip^TS 2 ) — > T^ip^TS 2 ) whose 
spectrum determines the stability properties of the harmonic map if> [23]. It is important to 
realize, however, that in H3.4[) is acting on Y, which is not a tangent section. So in (13.41) . 
is precisely the same operator defined above (|2.4p . but extended to act on sections of ^ -1 IR 3 . 
But is not self-adjoint (with respect to L 2 ) as an operator on t/* -1 !^, and self-adjointness 
of (the analogue of) is a crucial ingredient in Stuart's method. In fact 

J^Z = -AZ-(\if) x \ 2 + \if) y \ 2 )Z + A i ,Z (3.6) 

where the non-self-adjoint piece and its adjoint are 

A^Z = -2(lf> X ■ Z x +1Py- Zy)lP 

A\Z = -2{(iP-Z)Aifi + (iP-Z) x ifj x + (iP-Z) y iP y }. (3.7) 
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and we have adopted the analysts' convention for the Laplacian, that is, AZ = Z xx + Z yy . To 
remedy this deficiency, we make the following definition: 

Definition 3.1. Given a harmonic map ip : X — > S 2 , we define its improved Jacobi operator 
to be 

: r^" 1 !!) -> rxv" 1 !!), l^ = j^ + a\. 

Note that coincides with on T(ifj~ 1 TS 2 ), and hence maps tangent sections to tangent 
sections. Its principal part is the Laplacian, so it is elliptic, and it is manifestly self adjoint. 

Remark 3.2. Any section can be decomposed into tangent and normal components. As just 
observed, L$ maps a tangent section Z to the tangent section J^Z, so an alternative way of 
characterizing L$ is by specifying how it acts on normal sections, aip where a : X — > R. A 
short calculation, using harmonicity of if), yields 

L^(aif)) = -{Aa)ip - 4(a x ip x + a y ip y ). (3.8) 

It follows immediately that kerL^ = ker © (tp). Note that, in general, does not map 
normal sections to normal sections. 

Now, for any Y satisfying the pointwise constraint, 

A\Y = e 2 {\Y\ 2 Aij + 2(Y ■ Y x )^j x + 2(Y • Y y )ip y } =: j (3.9) 

and so, for any curve ip(t), if is a wave map then Y satisfies the PDE 

Y a + L^Y = k + ef (3.10) 

where j' = j + j. Note that j' has the same qualitative analytic properties as j (specifically, 
no higher than first derivatives of Y appear). 

We have yet to specify the curve ip(t) in M n . We do this by demanding that the error 
section Y should at all times be L 2 orthogonal to T^M n . In this case (S = T 2 , n > 2), 
T^M n = ker so, in terms of the local coordinate system q on M n provided by Proposition 
12. 1[ this amounts to requiring 

(Y,%) = 0, / u = l,2,...,4n. (3.11) 

We convert this into an evolution equation for q by differentiating the orthogonality constraint 
( 13.1 ip twice with respect to time and using ( I3.10p . 

(-L^Y + k + ef, ?/v> + 2eY t , ^ u q» + e 2 (Y, Vva)?V + e 2 (Y, Vv)<f = 0, (3.12) 

where an overdot denotes differentiation with respect to r. Now is self adjoint and £ 
ker C ker L^, so this equation simplifies to (k,ip^) = 0(e), or, more explicitly, 

<f + GlxW? = eh"(e, q, q, Y, Y t ) + e 2 ^(Y, ^ x }q x (3.13) 
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where 7 and G are the L 2 metric and its Christoffel symbol in the coordinate system q, (12. 7p . 
(12. 8p . and the function h is 

^ = 7^{(^,^A)g A + £(F,^ Ap )gV + ^((^A-^)y',^)gV 

+e 3 ((|F t | 2 -|F a; | 2 -|F,| 2 )y,^)}- (3.14) 

Taking the formal limit £ — > 0, ( I3.13P reduces to the geodesic equation on (M n ,7), as one 
would hope. 

To summarize, if is a wave map, and q(t) is a curve in M n such that Y = e~ 2 (4> ~ rfio)) 
satisfies the orthogonality constraint (13.1 ip at all times, then (Y, q) satisfies the coupled system 

Y tt + LY = k + ef, <f + G^q v q x = eW + e VC, $ uX )? (3.15) 

and the pointwise constraint (13. 3p . Conversely, if (Y,q) satisfies the constraints (13. 3p and 
(13. lip and the coupled system (13.151) . then = ip{q) + e 2 Y is a wave map. Our goal is to 
prove that (13.151) with fixed initial data q(0) = q , q(0) = q±, Y(0) = 0, Y t (0) = has solutions 
with ||y||c?i bounded uniformly in e for times of order It follows immediately that, in the 
limit £ — y 0, (j){r/e) converges uniformly to a curve ip(r) in M n . In the course of the proof, we 
will simultaneously show that ip(r) is the geodesic with initial data go^i- 



4 Analytic prelimaries 

In this section we set up the function spaces we will use, and collect some standard functional 
analytic results which we will appeal to repeatedly. More details can be found in [2], and 
references therein. Let H fc denote the set of real-valued functions on S whose partial derivatives 
up to order k are square integrable. This is a Hilbert space with respect to the inner product 

(f,g)k= Yl I D «fD a g (4.1) 

where a is a multi-index taking values from {x, y}, \a\ is its length and D a = d ai d a2 ■ ■ ■ d a ,,, 
so -D(a;,a;,j/) = d 2 d y , for example. We denote the corresponding norm by || • ||&, 

II/IB = </,/>*• (4-2) 

Let H h = H fc © H fc © H fc , the space of R 3 -valued functions on S whose components are in H fc . 
This is a Hilbert space with respect to the inner product 

(Y, Z) k = (Fi, + (Y 2 , Z 2 ) k + (Y 3 , Z 3 ) k (4.3) 

whose norm will again be denoted || • ||fc. We adopt the convention that || • || = || • || and 
(•) ') — (") ')o; that is, undecorated norms and inner products refer to L 2 . We will frequently, 
and withoutfurther comment, use the Cauchy-Schwarz inequality 

(Y,Z) k < \\Y\\ k \\Z\\ k (4.4) 
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and the trivial bound ||Y||fc > ||Y||fc' if k > k'. 

In the sequel, we will prove existence of a solution of the coupled system (13.151) with 
Y G H k , Y t G H , for k = 3. This choice of is motivated by the following fundamental 
fact about H fe on a compact 2-manifold. 

Proposition 4.1 (Algebra property of H fc , A; > 2). T/ie Banach space (H fc , || • ||&) zs a Banach 
algebra for all k > 2. T/iai is, if f aY\ k and g G H fe t/ien G H fc ; and t/iere exists a constant 
otk > 0, depending only on E and swca iaat H/gllfc < a fc||/||fc||fi , IU- 

It follows directly from this that, if (Y, 3^) G i? 3 © i? 2 , then the nonlinear term j' in the 
coupled system is in H 2 , and \\j'\\2 can be bounded by a polynomial in ||Y H3, \\Y t \\2 (the point 
being that j' contains no derivatives of Y higher than first, and Y x ,Y y , Y t are all in H 2 ). This 
is crucial, not only for proving the local existence result for ( 13.151) . but also in later sections 
where we prove that ||Y||3 is controlled by (LY, LLY), and make energy estimates for the 
solution. So the fact that we choose k = 3 is not just motivated by a desire to get strong 
bounds on the error section Y; the method will not work for any lower k. Indeed, to uniformly 
bound Y on S, it would suffice to control ||Y|| 2 , as we have the following Sobolev inequality. 

Proposition 4.2 (Sobolev inequalities). Let C k denote the Banach space of continuous maps 
S — t- M 3 with the usual norm \\Y\\ C k = sup{\D a Y(p)\ : \a\ < k,p G X}. Then H 2 C C° , 
H 3 C C 1 , and there is a constant a > 0, depending only on E, such that \\Y\\ C k < a||F||fc + 2 
for all Y G H 2 , k = 0, 1. (More briefly, the inclusions 1 : H 2 — > C° and 1 : H 3 — > C 1 are 
continuous.) 

In later sections we will need to bound ||L^F||fc in terms of ||F||fc + 2. Of course, since is 
a linear second order operator we have trivially, for all q G K, the upper bound 

\\L^Y\\ k < C\\Y\\ k+2 (4.5) 

where C is a constant depending only on E and K. For a lower bound, we use the fact that 
is elliptic. 

Proposition 4.3 (Standard elliptic estimate). Let D be an elliptic linear differential operator 
of order r acting on sections of a vector bundle V over E. Then there exist constants a^, fik 
depending only on E and k, such that 

||L>Y|| fc + a fc ||y|| > k \\y\\k+v 

If we consider only sections which are L 2 orthogonal to kerD, the same inequality holds with 
ctk = 0. 

The reason for quoting this result in the context of a general vector bundle V over E is 
that we will want to apply it to both the ordinary Laplacian on V = M?_, and the classical 
Jacobi operator on V = ip^TS 2 C IR 3 ,, where the H k norm is defined by inclusion. 
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5 Local existence theorem 

Theorem 5.1 (Local existence for the coupled system). Consider the coupled system $3.15]) 
with initial data q(0) = q G K, g t (0) = sq x G R 4n , Y(0) = Y G H 3 , Y t (0) = Y 1 G H 2 such 
that 



dist (go, dK) > d, 



|gi|,||r || 3 ,||Fi|| 2 <r 



where T,d are positive constants. Then there exist constants C(K) > and T(K,T,d) > 
such that for all e G (0,C(K)/\/r), this initial value problem has a unique solution on [0, T] 
with 

q G C 3 ([0,T},K) 

Y G C°([0,T],H 3 ) nC\[0,T],H 2 ) D C 2 ({Q,T], H 1 ). 

If the initial data are tangent to the L 2 orthogonality constraint 113.11]) and the pointwise 
constraint 113.3]) then the solution preserves these constraints. 

We will prove this using Picard's method: we iteratively define a sequence (q l , Y % ) G 
C°([0, T],K x H 3 ) which converges to a solution of the initial value problem. To establish that 
the iteration scheme is well-defined and convergent, the following standard energy estimate 
for the driven wave equation is key: 

Theorem 5.2 (Existence and energy estimate for the wave equation). The driven wave equa- 
tion on [0, T] x £ 

Y tt — AY = E 

with 5 : [0, T] x S — y M 3 smooth and smooth initial data Y = Y(0) and Y\ = Y t (0) has a 
unique global solution. The solution is smooth, and there exists an absolute constant c(S) > 1, 
depending only on the choice of torus E, such that 



max{||y t (t)|| 2 , \\Y(t)\\ 3 } < c(S)e* <^ ||Fi|| 2 + ||Fo|| 3 + 



n [s 



ydS 



Proof. Existence, uniqueness and smoothness follow from [7]. Let E(t) = \\Y t 
Then 



l + \\Y 



E'(t) = 2(Y t ,AY + E) + 2(Y,Y t ) 1 = 2(Y t ,S)o + 2(Y t ,Y) 

< 2||y t ||g + ||y||g + ||H||g<2£;(t) + ||r'" 11 - 



dt 



e- M E(t)) < e 
=> E(t) < e 



-2t||'-' 



o<l|S(*)||S 



E{0) + 



l^llo 



(5.1) 



(5.2) 



Now consider Z = AY. This is also smooth and satisfies the wave equation with source A: 
Applying the above estimate to Z yields 



|Ay t ||g+||AY||? < e 



2/ 



|Ay 1 ||2 + ||Ay ||? + 







/'l|AS||g 


< 2e 2t 


h 





mil + ll^c 



o ni + 



ii 



Since A is an elliptic operator, there exist positive constants o^, 0k depending only on k and 
£ such that 

||Ay||g + a fc ||y||g>^||y||2 +2 , 

by the standard elliptic estimate, Proposition 14.31 The result immediately follows. □ 

To prove Theorem 1 5. II we must first write the coupled system (13.151) as an explicit evolution 
system (note that both equations have q on the right hand side). Let X = M. 4n x M 4n x H 3 x H 2 
given the norm ||(g,p, Y, Z)\\x = max{|g|, ||y||3, ll^lh}- The e dependence of the norm 

is chosen so that ||(0, q t , 0, 0)\\x — <?• Given any T > let X r = {(q,p,Y, Z) G X : q G 
K, ||(0,p, y, Z)\\ < 8c(E)T} where c(S) > 1 is the absolute constant obtained from Theorem 
15.21 Note that Xr is a closed subset of a Banach space, and hence is a complete metric 
space with respect to the metric induced by || • Consider the matrix valued function 
M : R x I r -> End(M. 4n ), 

M (e, q, Y)* = 51 - eV A <y Vav>- (5-3) 

Since the matrix (7^ A ) is postive definite, K is compact, all g-dependence is smooth, and 
||y|| < 8c(S)r there exists a constant c(K) > such that M : [0,c{K)/VT] x X r ->■ 
GL(An, R) and M~ l : [0,c(K)/VT] x X r ->■ GL(4n, E) is C 1 and bounded. Hence, for all 
£ G [0, r)], where = c(K)/y/T the coupled system can be rewritten 

to = e 2 f{e,q,q t ,Y,Y t ) (5.4) 
y«-Ay = g(e,q,q t ,Y,Y t ) (5.5) 

where 

/(e,g,g t ,y,y t ) = M~ 1 (e, q, Y)(—G(q, q t , q t ) + eh{e, q, e~ l q u Y, Y t )) 

G(q,u,vr = G» x (qXv x (5.6) 

g(e,q,q t ,Y,Y t ) = -B^Y - ~ + ef(e, q, e~\,Y, Y t ), (5.7) 

B^, denotes the first and zeroth order piece of L^, so = —A + B^, explicitly 

b^y = -(i^,i 2 + i^i 2 )y-2(^-y, + ^-y/)^-2(^-y)A^ 

-2(V ■ y),^ - 2(^ • Y) y ip y , (5.8) 

and ft, and f = j + j are as defined in (I3.14p . (13. 5p . (13. 9p . 

It is convenient henceforth to consider e as a fixed parameter in [0, £*(K, T)] and supress 
the dependence of /, g on e. The proof of existence will use Picard's method, which requires 
that /, g be bounded and Lipschitz on X r . This follows quickly from the following proposition, 
whose proof is straightforward but lengthy, and so is deferred to the appendix: 

Proposition 5.3. The functions f,g are continuously differentiable maps f : Xr —> IR 4n and 
g : X r ->■ H 2 . Their differentials df : X r ->■ C(X,R 4n ), dg : X r -> C(X,H 2 ) are bounded, 
uniformly in e. That is, there exist constants Af(K,T), A g (K,T) > swc/i t/iat 

Id/^l < A/||a;||x, Hd/a^lh < A g |M|x 

for all x G X r , uj G X . 
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Note that, for Banach spaces B, C, C(B, C) denotes the space of bounded linear maps B 
C, which is itself a Banach space with respect to the norm H^U^^cr) = sup{||5 , (z) Hc/IMIb : 
x G B,x 0}. 

Corollary 5.4. The functions f : Xr — > M 4n and g : X-p — > H 2 are Lipschitz and bounded, 
uniformly in e. That is, there exist constants A/, A g , Cf, C g , depending only on K and Y , such 
that for all x, x' G Xr, 

\f(x)-f(x')\ < A f \\x-x'\\ x , \\g{x)-g{x')\\ 2 < A g \\x-x'\\ x , \f{x)\ < C f , \\g(x)\\ 2 < C g 

Proof. Let x\,x 2 G Xr. Since Xr is convex, the curve x(t) = x\ + t(x 2 — xi), < t < 1 
remains in Xp. Hence 



l/(*i)-/(* 2 )| 



/ df xit) (x 2 - xi) dt < / A f \\x 2 - XxWx dt = A/||a;i -x 2 ||x-(5.9) 
Jo Jo 



From the definition of / one sees that f(q, 0, 0, 0) = for all q G K. Hence, for all x = 
(g,p,y,Z)GX r 

|/(ar)| = |/(g,p,y,Z)-/(g,0,0,0)|<A / ||(0,p,y,Z)|U<8c(S)rA / . (5.10) 

The proof for g follows mutatis mutandis. □ 

To establish uniqueness of the solution, and to show that q is three times continuously 
differentiable, we will need the following extension property of df and dg, whose proof is also 
deferred to the appendix: 

Proposition 5.5. The differentials df : X r ->■ £(X,M 4n ) and dg : X r -»■ C(X,H 2 ) of f and 
g extend continuously to maps df ext : X r — > C(M. 4n x R 4n x H 1 x L 2 , M 4n ) and dg ext : X r — > 
C(M. 4n x M 4n x if 1 x L 2 ,L 2 ), bounded by Af and A g respectively. 

The rest of this section is devoted to the proof of Theorem 15.11 Let T > be chosen such 
that 

1 d r 

T<log2, T < T<— — — — — — , T< 



e ~ e(2T + C f {K,T)Y ~ eC f (K,Ty 
VT < , T< 7 , Vf < . (5.11) 

- c g (K,ry -4eA f (K,ry v - s c (s)a 9 1 ] 

Given a complete subset B of a Banach space with norm || • \\b, denote by CtB the space of 
continuous maps [0,T] — > B equipped with the sup norm ||6||| = sup{||6(t) ||b : t G [0, T]}. 
(CtB, HI • HI) is itself a complete subset of a Banach space. 

5.1 Definition of the iteration scheme 

We will produce a sequence uj % G CVX r converging to a solution of the initial value problem 
for the coupled system. Choose and fix 5 G (0, T/4), and let Yq, Y{ G C°°(E, M 3 ) be sequences 
such that 

\Y^-Yoh<L ||Y 1 i -yi|| 2 <^. (5.12) 
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Such sequences exist since C°° is dense in H k for all k > 0. Let u° = (go, £Qi, Yq, Y®), which 
is constant in t and smooth on E, and trivially lies in CtXy- Given u l , we define the next 
iterate to be the solution of the initial value problem g* +1 (0) = qo,q l t +1 (0) = £Qi,Y l+1 (0) = 
F m ,y/ +1 (0) = Y; +1 for 

e 2 f(tu*) 



& 1 



Y 



i+l 



AY 



i+l 



(5.13) 
(5.14) 



We must first check that the sequence uf is well defined. So, assume that u l is smooth and 
lies in CtXy- Then 

q l+ \t) =q + et gi + e 2 jf QT f(u\r))dr^ ds, 



(5.15) 



which exists since / o uj l is continuous. Now 

\q l+ \t) - g | < eT\ qi \ + ^e 2 T 2 C f (K,T) < eT (t + 
by our choice of T, so q l+1 (i) remains in K. Further, 



C f (K,T)\ < d 



e Iff* (*)l<l?i|+£ / \f(u l (s))\ds < T + TeC f (K,T) < 2T < 8c(S)r 



(5.16) 



(5.17) 



by our choice of T. Turning to Y l+1 , we see by inspection that if Y\ q l are smooth, then g(oj l ) 
is smooth, so the solution Y t+1 (t) exists, is unique and smooth, by Theorem I5.2[ which also 
yields the energy estimate 



max{||Y" t 



i+l 



\Y i+1 (t)\\s} < C^)e t {\\Yt 1 h+\\Y^% + V~tC 9 (K,T)} 

< 2C(E){2(T + 5) + VTC g (K,T)} <8C(E)T (5.18) 



by our choice of T and 5. Hence, if u l is smooth and in Cr-Xr, so is We have already 

observed that u° is smooth and in CrXr, so, by induction, the sequence u % G CtXy is well- 
defined. 



5.2 Convergence of the iteration scheme 

We will now show that u l is Cauchy, and hence converges in CtXt- From (I5.15P one has 



q i+1 (t)-q\t)\ 



by our choice of T. Similarly 



e- l \qt + \t)-<im\ = b 



t PS 



JO 



(f(u\r))-f(u*-\r)))drds 



< i-^A^r)!!^-^- 1 !!!^^^-^ 



i-l| 



{f{u\r))-f{u*-\r)))dr 



< eTAf(K, T)\\\lu 1 — u; || | < 



i-li 



(5.19) 



(5.20) 
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Now Z = Y l+1 — Y l satisfies the wave equation with source g(u) 1 ) — g(uo l *) and small smooth 
initial data ||Z(0)|| 3 , ||Z t (0)|| 2 < 5/2 i -\ Hence, by Theorem E2\ for each t G [0,T], 
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max{||Z(f)|| 8 JZ t (f)|| a } < c(E) C *^+(/ H^))-^*- 1 ^))!!^ 



i 

* \ 2 



< 2 C (E)|^j + VrA s (K,r)|||a; i -a; i - 1 | 

< ^ + l||| w *_ w *-i|||. (5.21) 



Assembling these inequalities, one sees that 



|o/ +1 -a/||| < Im^-^- 1 !!! + - (5.22) 



4' 



where a = c(£)5/8. It follows that 



I - wl I < ^III^-^III + ^t, (5.23) 

and hence, for all k > 1, 

|||u/+ fe -u/||| < V \\\u i+j -u^Wl < -|||o; 1 -a; |||y- + -y-. (5.24) 
in in — /_^\\\ in — ^ i n '"Z— ^47 2 l ^-^ 13 

3=1 j=l j=0 

Hence u l is Cauchy with respect to 1 1| • || |, so uj % — > uj = (q,p,Y,Z) G CxX r . 

5.3 The limit solves the initial value problem 

We have established that 

q i -> q in C T if, (5.25) 

gj 4 p in C T R 4n , (5.26) 

F' 4 F in C T H 3 , (5.27) 

->■ Z in C T H 2 . (5.28) 

Now, for all i, 



MO)-(go,egi,y ,yi)||x < lk(0) -u/(0)|| x + lk(0) - (g , egi,^,^ 



< |||u;-u/||| + -r ->■ (5.29) 
2* 

as z — > oo. Hence w(0) = (qo, eqi, Yq, Y\), that is, the limit has the correct initial data. 

We will now show that the limit solves the coupled system and has the differentiability 
properties claimed. Let Y(t) = Y + J^Z(s)ds. Note that Y is manifestly in C 1 ([0, T], H 2 ), 
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with derivative Y t = Z. Now 

ft 



\\\Y*-Y\\\ CtH 2 = \\\Y i (0)-Y o + ( {Y}(s)-Z{s))ds\\\c T m 

Jo 

< i+T\\\Yj-Z\\\ CTH ^0 (5.30) 

as % ->■ oo. Hence, Y i ->■ F in C r iJ 2 . But F< ->■ F in C T # 3 , hence also in C r # 2 , so F = F. 
Hence, F G ^([O, T], H 2 ) and F = Z. 

Consider Y^ k — Yfi.. This is smooth, and satisfies the wave equation with source g(u l+k ~ 1 ) — 
g(w l_1 ). Hence 

r t i +fe -^iii = iiA(F i+fc -F i ) + ( ? (^ +fc - i )-( ? (^- i )ii 1 

< 2\\Y^ k -Y% + \\g^ +k ' 1 )-g(co^ 1 )\\ 2 

< 2\\Y l+k - Y% + h g (K,T)\\u i+k - 1 - cu^Wx 

=> Ill^-^llleW < [2 + A g (K,r)]\\\u; i+k - 1 -u J i - 1 \\\c T x r . (5.31) 

Since w ! is Cauchy in C T X r , it follows that Yf t is Cauchy in C T H l . Hence Yf t — > W in CtH 1 . 
Let Z(t) = Fx + /J W. Note that Z is manifestly in ^([O, T],H r ) and = W. Now 

< 111^-^111^ + HI f{Y} t -W)\\\ CT m 

Jo 

< i+T\K-W\\\ CTm ^0 (5.32) 

as i ->■ oo. Hence F/ -» Z in CtH 1 . But F/ ->■ Z in C T H 2 , hence also in CtH 1 , so Z = Z. 
But F t = Z. Hence, F e C 2 ([0, T], tf 1 ) and F ti = IF. 

By similar reasoning, ^ ->m in CVIR 4n and q £ C 2 ([0, T], M 4n ) with qt—p and g it = m. 

We can now show that uj solves the coupled system: 

| \\Y U - AF - g(u>) || \ CT m < I - Yl\\ \ CtH , + | ||A(F - Y*) \\ \ CtII , + | \\g{u) - g^ 1 ) \\ \ CtH , 

< | ||F tt - F/j | CtH1 + 2| ||F - Y* 1| | Crif3 + | \\g(u>) - g^- 1 ) \\ \ CtH ,. 

(5.33) 

Now F4 ->■ F tt in CtH 1 , Y i ->■ F in C T H 3 , g : X r -> iJ 2 is continuous, and ->■ a; in C T X r , 
so — >■ in C T H 2 . Hence 

|||y tt -Ay-^(w)||| CTH i =0. (5.34) 

Similarly ||g tt - e 2 f{u) \\ Ct ^u = 0, that is, q tt = e 2 f{u). 

It remains to establish the higher differentiability of q. Differentiating the equation for 
tit 1 gives 

tit 1 = e'drfut. (5.35) 
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Now u l -> lo in CrXr, -> w t in CV(M 4ri x R 4n x H 2 x H 1 ), and <i/ ea; * is continuous, so 
g^ 1 ->■ say, in C T R 4n . Let m(t) = e 2 /(?o,£?i,*b,*i) + Jo 1 = ftt(O) + Note that 

m G C 1 ([0,T],M 4n ) and m t = & Then 

Mt-mWc^n = \\ [\qi tt -£)\\ OTR ±n<T\\qi tt -£\\ CTR . n ^0 (5.36) 

./o 

so q\ t m in C T R 4n . But -» g tt in C T R An , so g tt = m. Hence, g tt G C x ([0, T], M 4n ), as 
claimed. 



5.4 Uniqueness of the solution 

Assume that (q, Y) is another solution of (15.41) . (15.51) with the same initial data and regularity 
as (q, Y), and let (p, Z) = (q — q,Y — Y). Then (p, Z) satisfies the system 



Z tt -AZ=~(t), p u = e 2 T(t) 
with initial data Z(0) = Z t (0) = 0, p(0) = p t (0) = 0, where 

E(t) = g(u) - g@), T(t) = /(w)-/(a5) 
and w = (g, g t , F, w = (g, g t , F, Define 



= 11^)11? + ll^(Ollo + N 2 + -ib: 



'tl > 



(5.37) 



(5.38) 



(5.39) 



which, by the regularity properties of (p, Z), is continuously differentiable, and has E(0) = 0. 
Reprising the argument in (15.11) . which requires only that Z G H 2 and Z t G H , one sees that 



^(t)<2,E;(t) + ||S(t)||2 + |T(t)| 2 . 
Now, arguing as in the proof of Corollary 15.41 with oj(s) = u> + s(u — u) 



(5.40) 



o = \\9(.u) - j w o 



d^(s)(w - w) 



< A s max{|p|,e ^ptl, ||Z||i, ||Z t || } 

(5.41) 



by Proposition I5T51 Similarly |T(t)| <Aymax{|p|,£ 1 \pt\, \\Z\\i, \\Z t \\ }. Hence 

E'(t) < KE(t) 
where k = 2 + Aj + A 9 , whence it follows that 

d 



dt 



e~ Kt E(t) < 0. 



(5.42) 



(5.43) 



So e Kt E(t) is a nonincreasing, non-negative function which is zero at t = 0. Hence -E(t) = 
for all t, and we conclude that (p, Z) = (0, 0) for all t, that is, (g, F) = (g, Y). 
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5.5 Preservation of constraints 



Given the solution (q, Y) produced above, define for each \l G {1,2,..., 4n} 

a ^) = ( y i^>- ( 5 - 44 ) 

The L 2 orthogonality constraint is that a^(t) = for all /i, t. By construction, the coupled 
system implies that a M = 0. If the initial data are tangent to the constraint then a M (0) = a M (0), 
and hence a M (£) = for all t € [0, T]. 

Similarly, given the solution (q,Y) produced above, define \ : [0, T] x X — > R by 

X = r-^(g) + ^ 2 |r| 2 . (5.45) 

The pointwise constraint is that x = everywhere on E. Note that x(t) G H 3 for all t G [0, T] 
so x(^) : S — > R is continuous. Assume that x(0) = and x*(0) = 0, that is the initial data are 
tangent to the constraint. A straightforward, if lengthy, calculation using the coupled system 
and the harmonic map equation for ip shows that x satsifies the linear PDE 



Xtt - A X = 2e 2 | [2i> x ■Y x + 2iPyY y -\ip r \ 2 + Y-AiP- 2e^ T ■ Y t 

-e\\Y t \ 2 - \Y X \ 2 - \Y y \ 2 )] X + {Y ■ ^ x ) Xx + (Y ■ ^)x,} 
=: ax + biXx + hx y (5-46) 
where a(t) G H 2 , h(t),b 2 (t) G H 3 for all t. Let E(t) = \\x(t)\\j + ||Xt(t)Ho- Then 

E '(t) = 2(Xt,Ax + ax + b 1 Xx + b 2 Xy)o + 2(Xt,x)i = '2(Xt,X + aX + biXx + b2Xy)o 

< \\Xt\\l + \\x + ax + b 1 Xx + b 2 X y \\l < \\Xt\\l + 4x\\l < ^E(t) (5.47) 



where 



max 8(1 + ||a(t)|| 2 + ||6x(*)||^ + \\h{t)\\i). (5.4f 



Hence e~ Kt E{t) is a nonincreasing, non-negative function which is zero at t = 0, so E(t) = 0, 
whence ||x(^)l|i = f° r & U t- Since we already know that x{t) : S — )■ R is continuous, it follows 
that x — everywhere. This completes the proof of Theorem 15. II 



6 Near coercivity of the improved Hessian 

By repeatedly applying the local existence theorem, we can extend the solution of the coupled 
system whilever q remains in K and ||^||3 and ||Yt||2 remain bounded. So to prove long 

time existence, we must, among other things, bound the growth of ||^||3. The first step is 
to show that ||F||3 is controlled by the quadratic form (LY, LLY) or, more precisely, by the 
quadratic form Q 2 : H 3 — > R defined next. 
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Definition 6.1. For a fixed harmonic map i/j(q) we denote by Qi,Q 2 the quadratic forms 

Qf.H^R, Q 1 (Y) = J {|y x | 2 + |Fj 2 -(|^| 2 + |^| 2 )|F| 2 -4(^-y a; + ^-n)^-^} 
Q 2 :H 3 ^R, Q 2 {Y) = Q X {LY). 

Note that both Q\ and Q 2 are continuous, and that Qi{Y) = (Y,LY) for all Y G H 2 . It is 
also convenient to define the projection map P : H k — > H k 

P{Y) = Y — (ip ■ Y)ip 

which pointwise orthogonally projects Y{p) to T^ P )S 2 . 

Lemma 6.2. For all Y G H 2 , Q^Y) > Q 1 (P(Y)). 

Proof. Y = P(Y) + fij;, where / = ■ Y G H 2 . Now 

Q 1 (Y) = (Y,LY) = (P(Y),LP(Y)) + 2(f1>,LP(Y)) + (f1>,L(fil>)) 

= (P(Y),LP(Y)) + (W,L(fi>)) (6.1) 

since L is self-adjoint and maps tangent sections to tangent sections. But, as we saw in Remark 

L(/V) = "(A/)V - 4(/^ + fy^y), (6.2) 

so(^,L(/^)) = -(/,A/)>0. □ 

If our error section Y were a tangent section, the results of [5] would immediately imply 
that Qi is coercive, that is, Q%(Y) > c(g)||Y~|| 2 , orthogonal to ker J: 

Theorem 6.3 (Haskins-Speight, |5J). There exists a constant c(q) > 0, depending continuously 
on q, such that 

Q^Y) > c(q)\\Y\\l 
for all Y G H 1 satisfying ip ■ Y = 0, L 2 orthogonal to ker J^q). 

Unfortunately, ip ' Y ^ in our set-up, but is small (of order e 2 ). This means we can only 
establish the following "near coercivity" property for Qi. This will suffice for our purposes, 
however. 

Theorem 6.4 (Near coercivity of Q±). There exist constants c(q),c(q) > 0, depending con- 
tinuously on q, such that 

Q 1 (Y)>c(q)\\Y\\l-s 2 c(q)\\Y\\ 1 \\Y\\l 
for all Y G H 2 satisfying the pointwise constraint Ii3. 3\) . L 2 orthogonal to ker J^t q ) ■ 
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Proof. By Lemma Qi(Y) > (P{Y),LP{Y)) = (P{Y),JP(Y)) since L = J on tangent 
sections. Given any Z e ker J, (Z,P(Y)) = (Z,Y + ftp) = (Z,Y) = 0, since Z is pointwise 
orthogonal to ip and Y is L 2 orthogonal to ker J. Hence P(Y) is L 2 orthogonal to ker J, and 
so, by Theorem 16.31 there exists a constant c(q) > such that 

Q 1 {Y)>c{q)\\P{Y)f 1 . (6.3) 

Now, since Y satisfies (13. 3ft . 

\\P(Y)\\l = \\ Y + U 2 \Y\ 2 n\>\\Y\\l-E 2 (Y,\Y\ 2 i,) 1 (6.4) 
and, by the algebra property of H 2 (Proposition 14. ip . 

(Y,|Y| 2 ^)i < ||K||i|||K|Vll2<C||y|| 1 ||Y||2||V'||2 (6.5) 

where C > is a constant depending only on E. Combining (16. 3p . (16. 4p and (16. 5p . and noting 
that H^lta depends continuously (in fact smoothly) on q, the result immediately follows. □ 

Theorem 6.5 (Near coercivity of Q 2 )- There exist constants c(q),c(q) > 0, depending con- 
tinuously on q, such that 

Q2(Y)>c(q)\\Y\\l-e 2 c(q)(\\Y\\l + e 2 \\Y\\t) 

for all Y G H 3 satisfying the pointwise constraint Ii3. 3\) . L 2 orthogonal to ker J^t q ) ■ 

Proof. Recall that the pointwise constraint (13. 3 p is equivalent to \ip + e 2 Y\ = 1. The set of 
smooth maps S — > S 2 is dense in the Banach manifold of H 3 maps S — > S 2 , ip and ker J 
are smooth, and Q2 H 3 — > R is continuous, so it suffices to prove the inequality in the case 
that Y is smooth. So, let Y be smooth, L 2 orthogonal to ker J and satisfy the pointwise 
constraint (13. 3p . Define the smooth section Z = LY and the smooth real functions a = \Y\ 2 
and P = Z -Y. Since Y satisfies (I3.3p . and ip is harmonic, it follows that 

P = e 2 {Y ■ AY + \Y X \ 2 + \Y y \ 2 + 2(|^| 2 + |VvD|Y| 2 } • (6-6) 

In the following, Ci(q), 02(g), . . . denote positive functions depending continuously on q. We 
have the following elementary estimate, 

WMl < Ci(q)\\P\\l < e%(q){\\Y\\ 2 cl \\Y\\ 2 2 + \\Y \\ 2 c0 \\ Y\\ 2 3 } . (6.7) 
Applying the Sobolev inequalities (Proposition 14.21) gives 

||^||i<£ 2 c 3 (g)||Y||^ (6.8) 
We will also need to estimate [| 1 1 3 . Again, we have an elementary estimate 

H\i<c{\\Y\\ 2 cl \\Y\\ 2 2 +\\Y\\ 2 c0 \\Y\\ 2 3 } (6.9) 
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which, on appealing to Proposition 14.21 yields 



By Lemma 16.21 



Hl3<c||y||i 



Q 2 {Y) = Q 1 (Z) < Qi(P(Z)) = (P{Z), JP{Z)). 



(6.10) 



(6.11) 



It is in the last step that we have used the smoothness of Y (assuming only Y £ if 3 gives 
PiZ) £ H 1 , which is not sufficiently regular to make sense of JP(Z)). Since L is self-adjoint 
and ker J C ker L, Z = LY is automatically L 2 orthogonal to ker J, as is P(Z) = Z + ftip (since 
if> is pointwise orthogonal to anything in ker J). Hence, by Theorem 16.31 and the estimate ( 16. 8ft 



Qi(Y) 



> 



C4 



i^cMiwzwi-wzhimii} 



> c 4 (g){||Z||?-£ 2 C3 



Y\\l\\Z\\i} 



(6.12) 



We next estimate \\Z\\i = ||£y||i in terms of ||y||3. Note that Y is not L 2 orthogonal to 
ker L, since it has a component in the direction of if), so we cannot apply the standard elliptic 
estimate for L directly. We must decompose 

PiY) + (if) ■ Y)ifj = PiY) - \e 2 aifj, (6.13) 



y 



using (13.3p . and handle the two terms separately. Then, by Proposition 14.31 (for the lower 
bound on || JP(Y)\\x), and an elementary estimate (for the upper bound on || JP(Y)\\x), 



||JP(y)|| 1 --£ 2 ||L(a^||i 



< 



c 5 (g)||P(y)||3-£ 2 c 6 (g)|M|3< 

C5(g)||y|| 3 -£ 2 c 8 (g)||a||3< 
c 5 (g)||y|| 3 -e 2 c 9 (g)||y||l 



< 



l^lli 

\Z\\i 
\Z\\i 
\Z\\i 



< \\JP(Y)\\ 1 + -e 2 \\L(aiP)\ 
<c 7 (q)\\P(Y)h + e 2 c 6 



a 3 



<C 7 

<c 7 



Yh + e 2 c 8 
Y\\ 3 + e 2 c 9 



a\\3 

Y\\ 2 
1 II3 



(6.14) 



where we have used 
follows. 



.101) in the last line. Combining (16.121) and (16. 14(1 . the result immediately 

□ 



Remark 6.6. Since K is compact, we can replace c(q), c(q) in Theorems I6.4( [631 by global 
constants C, C > 0, under the extra assumption that q £ K. 



7 Energy estimates for the coupled system 

Having shown that Q 2 (Y) controls ||y|||, for small e, we must now bound the growth of Q 2 (Y) 
for a solution (q, Y) of the coupled system. We do this by establishing quasi-conservation of 
energies Ei, E 2 , related to Qi, Q 2 : 

Definition 7.1. Let (q, Y) : [0, T] — >■ K x H 3 with the regularity of Theorem 15. 1[ Associated 
to (q, Y) we define the energies Ei,E 2 : [0, T] — > R, 

Ei(t) = ImWl + \q x {Y), E 2 (t) = \\\{LY) t \\l + l -Q 2 {Y). 

Note that Ei is C 1 and E 2 is continuous. 
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Throughout this section we will use the following 

Convention 7.2. C will denote a positive constant depending (at most) on the choice of £ 
and K. c(ai,a 2 , ■ ■ ■ ,a p ) will denote a smooth positive bounding function of p non-negative 
real arguments, which may also depend (implicitly) on E, K, and which is increasing in each of 
its arguments. C , c will denote that the constant or bounding function depends, in addition, 
on the initial data qo, qi, Yq, Y\. The value of C, Co, c, Co may vary from line to line. 

Theorem 7.3 (Quasi- conservation of E\). Let (q,Y) : [0, T] — > K X H 3 be a solution of the 
coupled system with the initial data and regularity of Theorem \5.1[ Then 

E^t) < C + C(\q\ + \q\ 2 )\\Y(t)\\ + e [ c(\q\, \q\, \q\, \\Y\\ 3 , \\Y t \\ 2 ). 

Jo 

Proof. The solution satisfies (13.151) and has Y G H 3 , Y t e H 2 , Y tt e H 1 , so 

^ = (Y t , Y tt ) + (Y t , LY) + l -e(Y, L T Y) = (Y t , k + ef) + l -e(Y, L T Y) 
= ^-(Y,k}-e{Y,k T -h T Y) + e(Y t ,f) 
^E 1 (t) = E 1 (0)-(Y(0),k(0)) + (Y(t),k(t))+eJ S.(Y,^L T Y-k T ) + (Y t ,j')\ 

< Co + \\k(t)\\ \\Y(t)\\ + 2e [\WYWI + \\L T Y\\l + \\k T \\l + \\Y t \\l + . 

Jo 

(7.1) 

Now, only the first and zeroth order parts of L depend on time, so it is clear that 

IlL.Fllo^ClglllFll!. (7.2) 

Recall that k = —ip TT , so 

||^)llo<C||Mt)||c;o<C(|g| + |g| 2 ), (7.3) 

and 

\\k T (t)\\o < C\\k T (t)\\ c o < C(\q\ + \qf + \q\ 3 + |g| 2 ). (7.4) 

Finally, it follows immediately from the algebra property of H 2 (Proposition 14.11) that ||j'||o < 
lli'lb < c (|?|) ||^||3) ll^lh), and the result directly follows. □ 

We will need a similar result bounding the growth of E2(t). Formally, this is obtained by 
applying the argument above with Y replaced by LY (which formally solves a PDE of the 
form (LY)u + LiLY) = Lk + 0(e)). Unfortunately, this argument is not rigorous since Y is 
insufficiently regular to make sense of expressions like LY U (recall Y u is only H l ). 

Theorem 7.4 (Quasi-conservation of E 2 ). Let (q, Y) : [0, T] — > K X H 3 be a solution of the 
coupled system with the initial data and regularity of Theorem \5.1[ Then 

E 2 (t) < Co + C(\q\ + |<z| 2 )|mt)|| 2 + ^ / c(l?l. til Ma, ll^lh)- 
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Proof. By uniqueness, (q, Y) must arise as the limit of an iteratively defined sequence of 
smooth functions (q l , Y l ), as constructed in the proof of Theorem 15. II Recall that the sections 
Y i satisfy the PDEs 

Yjt +1 - AY l+l + B'Y* = Jfe* + e(j'Y (7.5) 

where B denotes the first and zeroth order piece of L (so L = — A + B) and the superscript 
i on B\ k\ (J 1 ) 1 denotes that the quantity is evaluated on the iterate (q\Y l ). Recall also 
(Y\ Yj, Yi t ) -> (Y, Y t , Y tt ) in C°(H 3 © H 2 © H l ). Now, for each i define 

Z i = -AY i+1 + B { Y { 
E\t) = \\\Zl\\ 2 + \{Z\-AZ l + B l Z l ). (7.6) 

Z Zi 

Each E* : [0,T] ->■ R is smooth, E i ->■ £ 2 uniformly on [0,T], Z* ->■ Z = LY in i/ 1 and 
Z t * -> (LF)t in L 2 . Note that Z* satisfies the PDE 

Zj, - AZ* + = B i (Z i - Z l ~ l ) + U + ef (7.7) 

where k = —Ak l + B % k % ~ 1 and f = —A(j') 1 + £> l (j') l_1 . Since B l is self-adjoint, one has 



1 



= <z;, b*(z* - z^) + fc 1 } + e {(zif ) + -<z\ 

= ^(Z i ,B i (Z i - Z 1 - 1 ) + jfe*) - (B^^Zl - Z\- x ) 

+e{(Z\ B^-tf - Z*- 1 ) + k\) + (Zlf)}. (7.8) 
Integrating this from to t and taking the limit i — > oo yields 

E 2 (t)-E 2 (0) < C + \\Z(t)\\ \\k(t)\\ + Ce I {||Z|| (||A; T ||o + |g|||^||i) + ||^||o||j||o} 

Jo 

(7.9) 

where we have used the facts that Z\ — Z' l f l — > in L 2 , Z % — Z l ~ l — > in H 1 , k l — > k = Lk in 
H 3 , k\ ->■ fc T in # 3 and j* -> j = Lf in L 2 (since L : H 2 -> L 2 and j : Kx R An xH 3 xH 2 ^ H 2 
are continuous). Now, we have the elementary estimates 

\\Z\\k = \\LY\\ k < C\\Y\\ k+ 2, k = 1,2, 
H^Ho = ||Ly t + eS T y|| <C7(||y t || 2 + e|g|||y|| 1 ) > 
\\k\\o < C(\q\ + \q\ 2 ), 

llillo < ^(|g| + l<?1 2 + l<?l 2 + |g| 3 ) 
lli'llo < c\\j'\\ 2 <c(\q\,\\Y\\ 3 ,\\Y t \\ 2 ). 

Combining these with (17.91) . the result follows. □ 
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8 Long time existence and proof of the main theorem 



Throughout this section we choose and fix q G K and q\ G M 4n , and denote by (q, Y) the 
solution of the coupled system (13. 15ft with initial data g(0) = go, ?(0) = qi, Y(0) = Y t (0) = 0. 
By Theorem 15.11 provided e < C(K)/ y\q~i\, this solution exists at least for time t G [0, To], 
where T depends on the initial data, but is independent of e. Moreover, the solution is unique, 
has the advertised regularity, satisfies the pointwise and L 2 orthogonality constraints (I3.3p . 
(13. lip , and obeys the energy estimates of section [71 Denote by g*(r) the geodesic in (M n ,7) 
with the same initial data, g*(0) = go, g*(0) = q±. Note that g*, considered as a function 
of rescaled time r, is independent of e. Since geodesic flow conserves speed 7(g*,g*), which 
uniformly bounds |g*| 2 on K, there exist r > 0, a > 0, depending only on the initial data, 
such that g* exists and has 

\q*\ < a , \q*\ < «o, dist (g(r), dK)) < d/2, (8.1) 

for all r G [0, r ], where d = dist (go, dK). Hence, the geodesic g* exists for time t G [0, £ _1 r ] 
which, for e small, exceeds Tq. Whilever g, g* both exist, we define s 2 q(t) to be the error 
between them, that is 

g = g* + e% (8.2) 

and 

M(s) = max {e 2 \q{t)\ 2 + W(t)\ 2 + |g"(t)| 2 + ||^(t)||§ + \\Y t (t)\\l} , (8.3) 

0<t<s 

where primes denote differentiation with respect to t. This function, which measures the 
total error in replacing the wave map cf> = ip(q) + e 2 Y with the geodesic ip(q*), is continuous, 
manifestly increasing, and has initial value M(0) = 0. Our next task is to bound its growth. 
Before doing so, we define another absolute constant (depending only on K and E), which 
will appear frequently in this section: 

a a = sup{|| 7 ^ A || : q G K, < fi, X < An}. (8.4) 

We will again use Convention 17.21 regarding bounding constants and functions. 

Theorem 8.1 (A priori bound). Whilever (q,Y) exists, and t < tq/e, and M(t) < e~ A a~ 2 , 

2 , . , ,2^2 , Co(M(f)) 



M(t) < C + C M(t)5 + ( £ 2 + £ t + eH 2 + eH 

'l-e 2 « a M(t)f 

Proof. We first derive the ODE satisfied by q. The curve g(r) satisfies the ODE 

q + G(q, g, g) = eh(e, g, g, Y, Y t ) + e 2 a(q, g, Y) (8.5) 
where G, h are defined in (I3.14p . (15. 7p . and 

a : K x R 4n x L 2 -> M 4n , a{q, v, Y f = 7 ^(^a, Y)v x , (8.6) 

which is smooth with respect to g and linear with respect to v and Y. The geodesic satisfies 
the ODE 

q* + G(g*,g*,g«) = (8.7) 
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with the same initial data. Substituting q = q* + e 2 q into (I8.5p . and using (18. 7p . we see that 
q satisfies 

q" = [G(g*, g*, g*) - G(g* + e 2 g, g*, g*)] + [G(g, g*, g* - G(g, g* + eg', g* + eg')] 

+£/i(e, g, g, Y, Y t ) + e 2 a(g, Y, g* + g"). (8.8) 

Now, g, by assumption, remains in K, and for all t e [0, T /e], g* remains in and |g*| < «o, 
so, since G(q,u,v) is Lipshitz with respect to g (on i^) and bilinear in (u,v), we have 

\q"\ < Ca 2 e 2 \q\ +C[e 2 \q'\ 2 + 600^1 +e\q\}e\h\+e 2 \a\ 

< C eM^ + Ce 2 M + e\h\ + e 2 \a\. (8.9) 

To estimate the h term, we note that it is smooth in g and polynomial in g and Y and its 
(first) derivatives, so in light of Proposition 14.11 we have (for e < 1) the crude bound 

\h(e,q,q,YY t )\ < c(\q\,\\Y\\ 3 ,\\Y t \\ 2 ). (8.10) 

Now, by the definition of M, 

\q\ < |g*| +e|?| < a + eM5, (8.11) 
so |/i| < Co(M). Turning to |a|, we have by linearity, 

|a(g,F,g; + g")| < a a ||il(|?.| + |g"l) < « a M^(a + M§) < c (M). (8.12) 

Hence, 

<ecb(M(*)). (8.13) 

Nowg^) = g-(0) + £g- = J^g 7 ', so 

W(t)\<£ [ c (M(s))ds <etc (M(t)) (8.14) 
Jo 

since Co and M are, by definition, increasing. Similarly 

e\q{t)\<£ I \q\<e 2 t 2 c (M(t)). (8.15) 
Jo 

We have now bounded the growth of all the q terms in M. To bound the growth of ||Y||3 
and || 1* || 2, we will use the energy estimates of section [7] and the near coercivity property of 
Q2 (Theorem 16.51) . But to do this, we need to control |g| and | g |, which appear in the energy 
estimates for E\(t) and E%(t), so we have not yet finished with the ODE for q. For q we have 
from ( I8.13P the obvious bounds 

|g| < |g'*| + |g"| < « + ec (M(t)). (8.16) 

For g we must work harder. So, for fixed e, let hi, i — 1, 2, 3, 4, denote the partial derivatives 
of h : K x R 4n x H 3 x H 2 — >■ M 4n with respect to each of its four entries. Similarly, let G 1: 
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a\ denote the derivatives of G, a with respect to their first entries. Differentiating (18. 5p with 
respect to r one finds 

? +Gi(q,q,q)q + 2G(q,q,q)q = ehiq + h 2 q + h 3 Y t + h 4 Y tt 

+e 2 ( ai (q, q, Y)q + a(q, q, Y)) + sa(q, q, Y t ) (8.17) 

where we have used the bilinearity properties of G and a. Inspecting the formula for h ( 13 . 14j) . 
we obtain, using (18. 161) and (18. lip a crude bound 

\Q\ < c (M) + co(M)\\Y tt \\+e 2 a a \q\\\Y\\ 

< c (M) + c (M)(\\LY\\ + \\k\\ + e\\j'\\) + e 2 a a M^\q\ 

< c (M) + c (M)(\\Y\\ 2 + \q\ + \q\ 2 + c(\q\, \\Y\\ 3 , \\Y t \\ 2 )) + e 2 a a M^\q\ 

< Co(M) +e 2 a a M^\q\. (8.18) 

Hence 

|g| < Co(M) , (8.19) 

1 - e 2 a a M2 

whilever M(t) < e~ A a~ 2 . We can now turn to bounding Y . 
Since Q 2 is nearly coercive (Theorem 16.51) . 

\\Y\\l<C{Q 2 (Y)+s 2 c(\\Y\\ 3 )}, (8.20) 

and, by the standard elliptic estimate for L (Proposition 14. 3[) 

\\Y t \\l < c{||Lri|| 2 + ||rj 2 }<c'{||(Lr) t || 2 + £ 2 ||L T r|| 2 + ||r t || 2 } 

< C{\\(LY) t \\ 2 + e 2 \q\ 2 \\Y\\l + \\Y t \\ 2 } (8.21) 
since the principal part of L does not depend on time. Adding (I8.20p and (18.211) . one sees that 

\\y Ha + M < C{^i(*) + + eMails) + e 2 |?| 2 ||*1i} < C{S x (t) + E 2 (t) + e 2 c (M)}, 

(8.22) 

where E\, E 2 are as in Definition 17. II Then, by Theorem 17.31 and 17.41 and the estimates ( 18. lip . 
dHZED, (Ell), 

M + ra < Co + (Co + ec (M))M^ +e f U+ C °[ M) x ) co(M) + e 2 c (M) 

Jo V l — e 2 a a M^/ 

< C + C M^ + ec (M) + £tCo(M) r , (8.23) 

1 — e 2 a a M2 

provided e < 1 and M(t) < £- 4 a" 2 . Combining (IQgj) . <KT^ . (EH} and ( Jgjgj) gives the a 
priori bound claimed. □ 

Theorem 8.2 (Long time existence). There exist > and r* > 0, depending only on the 
initial data, such that for all e G (0, £*) t/ie solution (q, Y) persists for all t e [0, r*/e] and aas 
M(t) bounded, independent of e. 
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Proof. Choose and fix a constant M* > 4«q so large that M* > Cq + CqM* , where Co is the 
specific constant, depending only on initial data, appearing in the a priori bound. Assume 
e > is so small that 



-2.1 _ d 1 ^ C(AT) 



^ < r, e < r, e < -, e < (8.24) 

n„.\/7 4M? " .U. T 



where, as before, d = dist(g ,9i^) and C(K) > is the constant quoted in Theorem 15.11 
Then, whilever M(t) < M„ ||r|| 3 < m), ||y t || 2 < M"|, 



and 



|£| < a + eM* < Mi, (8.25) 

e 2 \\q\\ <eM5 < £ M| < ~. (8.26) 
Hence, whilever M(t) < M*, the value of the solution (q,q t ,Y,Y t )(t) satisfies the conditions 



of the initial data for the local existence theorem 15. 1[ with T = M* and dist (q(t),dK) < |. 
Given the last condition on e, (|8.24|) . it follows that we may apply Theorem 15.11 and extend 
the solution for a time 5T > depending only on M* and d, independent of e. It follows 
that the solution persists for as long as M(t) < M*. Furthermore, by the first condition on e, 
f )8.24p . whilever M(t) < M* the solution obeys the a priori bound, Theorem 18.11 

For each e > let t £ = sup{t : M(t) < M*}. We claim that there exists > such that 
et £ is bounded away from zero on [0,5*]. Note that this immediately implies the statement 
in the theorem since then there exists r* > such that t £ > r^/e for all e G (0,5*), and the 
solutions exists, with M(t) < M* on [0, r*/e]. Assume, towards a contradiction, that no such 
exists. Then there is a positive sequence £j — >• such that e{t Ei — > 0. But then the a priori 
bound at time t £i gives (recall M is continuous, so M{t £i ) = M*), in the limit i — > oo, 

M*<C + C m|, (8.27) 
a contradiction, by our choice of M*. □ 



Proof of Main Theorem. By Theorem I8.2[ for all e G (0, £*] the solution exists for t G [0, r*/e], 
and coincides with if}(q*(r) + e 2 q{t)) + e 2 Y(t), with e|g(t)|, ||F(t)|| 3 , and hence (by Proposition 
4.2p ||y(t)||cfo uniformly bounded in t and e. The rescaled solution is 

£ : [0, rJxS4S 2 C M 3 , £ (r, x, y) = ^*(r) + £ 2 g(r/e), x, y) + £ 2 F(r/e, x, y), (8.28) 

and the geodesic with the same initial data is 

i:[0,r»]xS, ifj*(r,x,y) = ip(q*(T),x,y). (8.29) 

Now ip : K x X — > S 2 is smooth, hence uniformly continuous (since X x £ is compact). Hence, 
as e — > 0, £ converges uniformly on [0, r*] x E to f,. Furthermore, 

&l(T,x,y) = i) x (q*(T)+e 2 q(T/e),x,y)+e 2 Y x 
</>l(T,x,y) = i> y (q*(T) + e 2 q(r/e),x,y) + e 2 Y y 

<j> £ T {r,x,y) = (qi; + eq'(T/e))4>,(q*(T)+6 2 q(T/e),x,y) + eY t (8.30) 
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and |(f(£)|, ||Yj B (£)||c'o, ||l^(t)|| c o, ||Yt(£)||c° are bounded uniformly in t and e (again using 
Proposition 14.2)) . so (p £ x , (jf y , (p £ T converge uniformly on [0, t*] x E to ip* x ,iJ*y,4>*T- Hence, £ 
converges to in C 1 . □ 



Appendix: Analytic properties of the nonlinear terms 

The proof of the local existence theorem 15.11 makes fundamental use of certain basic analytic 
properties (Propositions 15.31 and 15.5)) of the right hand sides /, g of the evolution system 
f l5.4p .( l5~3]) . These properties are established by a long chain of elementary arguments which 
we sketch in this appendix. We regard £ as a fixed parameter in (0, 1) and choose a local 
parametrization ip : U x E — > S 2 of M n as given by Proposition 12.11 and a compact convex 
neighbourhood K C U . We begin by showing that the associated maps ^ : U — > H k , 
q t- y i/)(q, •) are smooth for all fcGN. 

Lemma A.l. Let k e N and / : [/ x E — >• E be smooth. Then F : U ->■ H fc , F(g) = /(<?, ■), zs 
smoot/i. 

Proof. As usual, we will denote partial derivatives with respect to q^ by a subscript /i. It 
suffices to show that F is everywhere differentiable, since all partial derivatives / MlA t 2 -/v are > 
like /, smooth maps U x E — > R. By the mean value theorem there exist t±, . . . ,t r G [0, 1] 
such that 

471 

\\F(q + p)-F(q)-f tM (q,.) P ql < \p\ 2 [ {(U(q + hp) - f,(q)f 

+(f^(q + t 2 p)-Ux(q)Y + --- 

• ■ ■ + (Uy-yiQ + t r p) - f„{q)f) (A.l) 

since / and all its partial derivatives are C 1 functions of q. But f^, . . . , f m - y are uniformly 
continuous on B$(q) x E, for 5 > sufficiently small, so 

limJ-\\F(q + p)-F(q)-^- p»\\ k = (A.2) 
p^o |p| dqf q 

as was to be proved. □ 

Now, by the definition of H k , ip : U — >■ if fc is differentiable if and only if ^ : f/ — > H fc is 
differentiable for i = 1,2, 3, so we immediately obtain: 

Corollary A.2. ^ : U — » g >->■ ^(g, •), smooth for all k 6 N. 

The error terms are j'(?> OS ^> ^t) an d ^(?, <L Y, Yt) where 

f(q,p,Y,Z) = 2p^^-Z)ij + e(\Z\ 2 -\Y x \ 2 -\Y y \ 2 W + ep^(^-tP u )Y 

-2e(^ ■ y x + ^ • y,)y + ^ 2 {|y| 2 a^ + 2(y • ya^ + 2(y ■ yjvj 

h(q,p,Y,Zr = Y u {(Z,^ uX )p x + e(Y,^ Xl/p )pY + e^ x -^ p Y,^}pY 
-2e{(ip x -Y x + ip y - Y y )Y, i, v ) + 2e 2 ((^ x ■ Z)Y, ^ u )p x 
+e 3 ((\Z\ 2 -\Y x \ 2 -\Y y \ 2 )Y,^)}. (A.3) 
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We will also need to consider the quantities 

M» v (q,Y) = 5^-e 2 Y X (q)(Y,^x u ) (A.4) 
A(q,Y) = -(\^ x \ 2 + \^ y \ 2 )Y-2^ x -Y x + i/j y -Y y )^-2^-Y)A^ 

-2(V> ■ Y) x ip x - 2(^ ■ Y) y ip y . (A.5) 

Let B denote the Banach space M. 4n x M 4 ™ x H 3 x H 2 with norm 

\\(q,p,Y, Z)\\ B = max{|g|, \p\, \\Y\\ 3 , \\Z\\ 2 }, (A.6) 

Bjj C -B denote the open set on which q G U and for each r > 0, By C denote the closed 
convex subset on which q <E K and ||(0,p, Y, Z)\\ B < Y. 

Proposition A. 3. The quantities defined above are smooth maps j 1 '. By — y H 2 , h ; By — y 

R 4n , M» v : By -> M, A: By -> H 2 

Proof. That j', /i define maps By — > H 2 and By — > M. An follows immediately from the algebra 
property of H 2 (Proposition |4J]). Now j' is a linear combination of terms formed by composing 
the maps 



dty k :Ux R 4n ->• H\ (q,p) 4 V^^M 



(A.7) 



which are smooth by Corollary IA.21 and the manifestly smooth maps 



H fc 4 fc , f 4 /e, 
H 2 xH 2 4H 2 , (f,g)^fg, (A.8) 



where ei = (1, 0, 0), e 2 = (0, 1, 0) and e 3 = (0, 0, 1). Hence j' is smooth. 

The map h is handled similarly, after noting that the inverse metric coefficients 7^ are 
smooth U — > M, the higher derivatives 

d 2 ^ fc : U x M 4n x M 4n 4 (g,Pi,p 2 ) 4 pfr&Mg, •) 

d 3 ^ fc :(/xl 4 "x M 4n x M 4n 4 # fc , (q,Pi,P2,p 3 ) ^ P»3?/W<7, •) (A.9) 

are smooth by Lemma IA.2I and, in addition to the maps in (|A.8[) above, the map 

#° x F° 4 R, (7, Z) 4 (7, Z> (A. 10) 

is manifestly smooth. 

That M^ v and A define smooth maps on By is clear, since their g-dependence is smooth, 
and they depend linearly on 7 (and are independent of p and Z). □ 
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We can now assemble these pieces to show that / and g are smooth functions on Xr, the 
space defined in section [5j To do so, we note that f = f°i and g = got where i : X r — > -Bs c (s)r 
is the linear isometry 

t:(q,p,Y,Z)^(q,s- 1 p,Y,Z) 

and f-.Bjv-t R 4n ,g : B r > H 2 are 

f(q,p,Y,Z) = M- 1 (q,Y)(-e 2 G(q, P ,p)+eh(q,p,Y,Z)) (All) 
g(q,p,Y,Z) = -A(q,Y)-^r(q,p,Y,Z)-^ u p^ + ej'(q,p,Y,Z), (A.12) 

and G is defined in (15. 6ft . The point is that f,g are defined as functions of q t (and q,Y,Y t ) 
on a space (Xr) with e-dependent norm, but it is more covenient here to think of them as 
functions of q = e~ l qt, on a space (-Br', T' = 8c(S)T) with fixed norm. Since i is a linear 
isometry, /, g are smooth, bounded, Lipschitz, etc. if and only if /, 7j are. 

Proposition A. 4. There exists = 0(1/ y/T) such that for alls G (0, £*), / : Br — >■ M 4n and 
~g : Br — > H 2 are smooth. 

Proof. Since 7 AtA (g), 4>\v(q) are smooth and X is compact, there exists > such that the 
matrix M(q,Y) is uniformly invertible on Br for all e G (0, £*). The components of M _1 
are rational in M^ v with denominator det M, which is bounded away from on Br- Hence, 
M _1 is smooth on Br, and the proposition follows immediately from the Leibniz rule and 
Proposition IA.3I □ 

Since / and 'g are smooth, they are certainly continuously differentiable. To complete the 
proof of Proposition 15. 3[ it remains to show that their differentials are bounded. 

Proposition A. 5. For all e G (0, £*), the derivatives d<7 : Br — > C(B,H 2 ) and df : Br — > 
C(B,M. 4n ) are bounded, independent of e. 

Proof. This is established by estimating the operator norm of the derivatives termwise. For 
example, the first term of j', 

J(q,p,Y,Z) = 2p^^-Z)^ (A.13) 

has derivative 

dJ (q , P ,Y,z) ■ {lp,Y : Z)^2q v p»[{^-Z)^ + {^^Z)^ (A.14) 

and so, by the algebra property of H 2 , for all (q,p, Y, Z) G -Br, 

\\dJ {qjPjY , z) (q,p,Y,Z)\\ 2 < C(K){\qM\\Z\\ 2 + \p\\\Z\\ 2 +\\Z\\ 2 \p\} 

< C(K){2T + T 2 }\\(q,p,Y,Z)\\ B (A.15) 

where C(K) > is a constant depending only on K. Hence, for all (q,p, Y, Z) G -Br, 

\\dJ {p , gX z)\\c(B,Hi) < C(K){2T + T 2 }. (A.16) 

The other terms of j', h and A are handled similarly. To bound dM _1 we bound dM and appeal 
to the Leibniz rule and uniform invertibility of M. Boundedness of the differentials of functions 
depending only on q and p is immediate by compactness and finiteness of dimension. □ 
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Finally, we turn to Proposition 15. 5[ which is equivalent to: 
Proposition A. 6. The differentials of the maps f : -Br R 4ra g : Br — > H 2 extend to maps 

d Jext . Br _^ £( R 4n x R 4n x jjl x j2 ^ R 4n) md d ~ezt g . ^ _^ £( R 4n x R 4n x #1 x L 2 , L 2 ) 

bounded by Af, A g respectively. 

Proof. This follows from explicit termwise computation. For example, the last term of h is 

m^q, Y, Z) = {(\Z\ 2 - \Y X \ 2 - |ig 2 )F, ^) 

whose differential at (q, Y, Z) is the linear map M 4n x H 3 x H 2 — > M 

dm M : (g, K, Z) ^ ((|Z| 2 - |F X .| 2 - \Y y \ 2 )Y, + ((\Z\ 2 - \Y X \ 2 - \Y y \ 2 )^, Y) 

-2((Y • ^)Y„ t) - 2((Y ■ ^)Y y , Y y ) + 2((y • Z) 

which clearly extends to a bounded linear map M 4n x iJ 1 x L 2 — > R. □ 
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